Geodesic and Path Motion in the Nonsymmetric Gravitational Theory 



in 



> 

in 



in 

0\ 



X 



J. Legare and J. W. Moffat 
Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 1A7 

(September 21, 1995) 

We study the problem of test-particle motion in the Nonsymmetric Gravitational Theory (NGT) 
assuming the four-velocity of the particle is parallel-transported along the trajectory. The predicted 
motion is studied on a static, spherically symmetric background field, with particular attention paid 
to radial and circular motions. Interestingly, it is found that the proper time taken to travel between 
any two non-zero radial positions is finite. It is also found that circular orbits can be supported at 
lower radii than in General Relativity for certain forms of motion. 

We present three interactions which could be used as alternate methods for coupling a test- 
particle to the antisymmetric components of the NGT field. One of these takes the form of a 
Yukawa force in the weak-field limit of a static, spherically symmetric field, which could lead to 
interesting phenomenology. 



C/3 ' I. INTRODUCTION 



The problem of particle motion in gravitational theory is a long-standing and complicated one. Yet, no problem 
could be more important in verifying the physical validity of a theory, since we measure a gravitational field by 
observing its effect on particle motion. 

Several studies have been written on the problem of motion in General Relativity (GR), both for test-particles and 
(f-) | for massive bodies. Einstein Q originally postulated that the four-velocity of a test-particle was parallel-transported 
along itself during the course of the particle's motion. He also showed that such a motion described a path of 
s * ' extremal length. In 1927, Einstein and Grommer showed that the field equations of GR did not allow for arbitrary 
motions of mass concentrations, rendering Einstein's original postulate of geodesic motion for test-particles spurious. 
Beginning in 1938, several workers sought to derive the equation of particle motion from the field equations, either by 
considering how the field equations governed the motion of singularities in the gravitational field Q , or by investigating 
O the matter-response equation, V ' V T^ = where is the stress-energy tensor, and its effect on the motion of mass 
concentrations ||,|| . In || , Wheeler discusses the potential problems associated with attempting to derive the equation 
5h i of motion from the gravitational field equations. In particular, he points out that describing concentrations of mass- 
if)' energy by singularities in the metric can render a theory incomplete. Meanwhile, approaching the problem from the 
point of view of the matter-response equation is troublesome, as different stress tensors yield fundamentally different 
equations of motion. As Wheeler points out, however, both of these methods are valid in the weak-field regime. 

The Nonsymmetric Gravitational Theory (NGT) involves a nonsymmetric, second-rank tensor field g^ v , interpreted 
as the gravitational potential. This tensor is refered to as the fundamental tensor. Although it is tempting to proceed 
in analogy with GR and associate with NGT a metric space, it is not clear a priori how this structure is to be identified. 
In old NGT (mathematically equivalent to Einstein-Straus Unified Field Theory, but with a different interpretation), 
it was found by means of a Cauchy analysis of the field equations (see and HH ; note that this work was performed 
in the framework of Einstein-Straus Unified Field Theory) that the theory contained three possibly different sets of 
characteristic surfaces, corresponding to three different metric structures (see (30.1)-(30.3) in It is not clear 

which should be chosen. 

In a recent version of NGT, also known as massive or finite-ranged NGT, the field equations have been modified so 
as to correct asymptotic consistency problems (see JTo| [l3| ) . In the linearized theory, this corresponds to endowing the 
antisymmetric sector of the gravitational field with a mass. Certain modifications to the compatibility condition have 
also been made (see, for instance, (10) of p2|) . These changes are sufficiently significant to cast doubt on the direct 
applicability of the precise results of and |8j|J. It is evident, however, that the metric structure of a generic NGT 
spacetime still would not be unique. Current research is exploring the possibility of generalizing the work described 

in § and®. 

Ignoring for the time being these unanswered questions, we will assume the metric structure is given by s^ v = g(a V ); 
the generalization of most of our results to some other metric structure amounts essentially to redefining certain 
metric functions. We will define g^ v by g^g^a = g vtl 9a^ = 5%, and s^ iV by s^s^ = 8%. Note that in the first of 
these, the order of the indices is important. We will also define a p „ = g^v]- A boldface character will be used to 
denote a tensor density: — ^J—gg^w 
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The problem of particle motion is no simpler in NGT than in GR. Despite its considerable historical value, an analysis 
of the type performed by Einstein et al. in S is not desirable due to the interpretational difficulties associated with 
introducing singularities into a field theory. On the other hand, the method of Papapetrou demands a priori 
knowledge of the stress tensor , or at the very least a physical interpretation of its components. In GR, the stress 
tensor is a symmetric tensor; this allows a simple physical picture of its components as being energy and momentum 
densities. In NGT, however, the stress tensor need not be symmetric. Although its symmetric components can 
be interpreted as in GR, is it unclear how the antisymmetric components should be interpreted. Furthermore, a 
discussion analogous to the debate on the uniqueness of the metric structure outlined above also applies to the stress 
tensor: it is not clear which form of the stress tensor (T^ v , T^ v or some other combination) should be interpreted 
as representing the energy and momentum densities associated with the gravitational field. It will nevertheless be 
informative to investigate the predictions of the matter-response equation for the unambiguous case of a structureless 
test- particle; this will be done at the end of Section HI. The methods of both Einstein || and Papapetrou [^j5| have 
been, to one extent or another, generalized to the case of NGT [ pT[ . 

Notwithstanding the criticisms of Wheeler, we will follow Einstein and simply postulate that the four-velocity of a 
test-particle is parallel-transported along itself during its motion. Section || discusses some generalities about particle 
motion and parallel-transport; we also introduce the notation used in the remainder of this paper. In Section III, 



we derive the first integrals arising from the equations of motion due to symmetries in the background spacetime. 
Section [V applies the results of the previous sections to the so-called geodesic and path equations for the special 
case of a static, spherically symmetric background spacetime. Radial and circular motions in a static, spherically 
symmetric field are discussed in Sections |V| and VI, respectively. In Section VII, we briefly demonstrate that the 
geodesic and path equations predict similar results in the weak-field regime; we derive the first-order correction to 
the geodesic equat ion c oming from the path equation, showing that this correction must be consistently ignored. 
Finally, in Section VIII, we discuss the possibility of modifying geodesic motion so as to include other couplings to 
the antisymmetric field. 



II. GEODESIC AND PATH EQUATIONS 

We define a connection D by its effect on a set of basis vectors e M : Dxe^ = CPX v ep, where C@ v are the coefficients 
of the connection, which need not be symmetric on [iv. We add to this the requirement that obey the Leibnitz 
rule, as well as commuting with contractions. By decomposing vectors onto this basis, we deduce that the effect of 
the connection D on vectors is given in component form by: D^v x — d^v x + C^v 9 . This establishes the operational 
use of D on vectors; the generalization to covectors and tensors of higher-rank is immediate. 

A vector u is said to be parallel-transported with respect to a given curve x(r) if its components at a point P' on 
the curve are proportional to its components at another point P on the curve, an infinitesimal distance away (see [ p5| , 
p. 13): D±u = <P(t)u, where x — dx/dr and 4>{t) is an arbitrary function of r. Here, r is some parameter for the 
curve x(r). In particular, if we take u = x, we may write in component form: u v D v vP = §{j)vP . The term on the 
right-hand side can be eliminated by affinely reparametrizing the curve, leading to the generalized path equation (see 
also 

u v D vU P = — + C^u v = (1) 

for some affine parameter r. Being nothing more than the definition of parallel-transport of a vector with respect to 
a given connection D, the generalized path equation is fundamental in itself; it is not clear that it could be derived 
from any more fundamental principle, say, a variational principle. The relation of the generalized path equation to 
paths of extremal length will be discussed in more detail below. 

We will take u 13 to represent the components of the particle four- velocity, and we will call ([!]) the equation of particle 
motion. 

Our definition of the generalized path equation is essentially that of Eisenhart (see p. 12). However, alternate 
but equivalent definitions of the generalized path equation exist in the literature. For instance, Lichnerowicz defines 
a path x(t) of the connection with coefficients C@ by: 

+ u« - + c»yv) u' = o, (2) 

(see [Q, pp. 248-249; there, the generalized path equation is referred to as a geodesic equation, regardless of the 
coefficients C@), where u — dx/dr is the tangent vector to the path. We can recover this definition from the 
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generalized path equation by multiplying ([!]) by u a and antisymmetrizing on a/3. Tonnelat (see |T^], pp. 145-148) 
and Hlavaty (see @, pp. 64-66) propose similar definitions. Although these are equivalent, (§) has the advantage 
of being independent of the parametrization of the curve x(t), and holds even for a non-affine parameter. For our 
purposes, the added benefit of being able to use a general parameter is not worth the increase in complexity associated 
with using a tensor equation of motion over a vector equation of motion. We will therefore use (Q) instead of (^J). 

Before proceeding any further, we must specify the coefficients C@. Although a multitude of connections can be 
formed in NGT, there are three that arise more commonly. The first is the symmetric Levi-Civita connection, whose 
coefficients are the usual Christoffel symbols: 

O = \^\d v s, x + - d xSflu ). (3) 

Being compatible with s M „, it is an appealing candidate for describing test-particle motion, as it leaves the magnitude 
of vectors, and hence the four-velocity, unchanged under parallel transport. The second is the unconstrained W- 
connection whose torsion vector appears in the NGT action pl|-p^| . The third is the T-connection, whose torsion 
vector is constrained to vanish. These last two connections are not independent, but are related by 

r^ = we u -hew v , (4) 



which ensures that the torsion vector r„ = 0. It can be shown (see 15 , pp. 12-13 and pp. 30-31) that the relation 
([|) is sufficient to guarantee that ([!]) will predict the same motion, regardless of which of these two connections is 
used. This can be seen immediately by introducing (|]) into (^). The problem is therefore reduced to studying motion 
predicted by a Levi-Civita connection and that predicted by a T-connection. 

Note that for a non-Levi-Civita connection, the solutions of ([!]) are not paths of extremal length [|l9| . By contracting 
(tih with s a /3 and collecting terms, we arrive at 



1 / d dL dL 



= ({£}-<&) , (5) 



2 \dr du a dx a 

where we have introduced L = s^u^u" as the magnitude of the vector u. We recognize the left-hand side of (||) as 
the Euler-Lagrange equation for the Lagrangian L, the very Lagrangian used when extremizing the length of a curve 
whose tangent vector is u. We conclude that if the right-hand side of (|J) vanishes, which is by no means assured, 
then the length of the path described by ([!]) will satisfy the Euler-Lagrange equation and the path will have extremal 
length. This occurs trivially when the connection appearing in (|l|) is the Levi-Civita connection. Other connections 
must be examined on a case-by-case basis. 

To be specific, we will refer to the case of (|l]) with the Levi-Civita connection as the geodesic equation; when using 
the T-connection, we will refer to this as the path equation. 



III. CONSERVATION LAWS 



We will now show that the generalized path equation contains first integrals of the motion that reflect a symmetry 
in the background geometry, allowing us to simplify the description of the motion. Interestingly, we will see that the 
first integrals are not always explicitly guaranteed to exist, leading us to question whether a conserved quantity is 
necessarily associated with every symmetry in the background geometry. 

Consider Xpu^ 3 , where Xp = sp a X a and where satisfies ([l]); X 13 are the components of an arbitrary vector. 
Take 



= u fi u v D u Xp + X p u v D v u f3 = u f3 u v D {u X p) - 



by direct calculation, 



D^Xrf = -£ X [*U + X P ({£} - C% 



where £x[$]fiv = X a d a s liv + Sa^d^X 01 + s^ a d u X a is the /x^-component of the Lie derivative of the symmetric part 
of the fundamental tensor along a flow generated by X. Therefore, 

= \£x[sUu»u» + Xpu»u" ({£} " <#,) ■ (6) 
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If X is a Killing vector, then £x [s] — and (||) reduces to 

^-W'(U}-C»). (7) 

If the connection coefficients C@ are taken to be Christoffel symbols, we see that (^) is sufficient to guarantee the 
existence of a conserved quantity. Two notable examples are X — d/dt and X = d/d<j), both Killing vectors of a 
static, spherically symmetric system. Writing s tt = 7 and = /3sin 9, we conclude from (Q) that Stpu^ =ji = E 
and s^pu^ = [3(f) sin 2 = J are conserved quantities, corresponding physically to the conservation of energy per unit 
mass and angular momentum per unit mass. 

The situation is radically different when the C@ are taken to be some connection coefficients other than the 
Christoffel symbols. In particular, take C^„ = T% u . In this case, we see from (^) that the fact that X may be a 
Killing vector is not sufficient to guarantee the existence of a conserved momentum; we must further require that the 
relevant connection coefficients be equal to the corresponding Christoffel symbols. A case in point is the Killing vector 
X = d/dt of a static, spherically symmetric system. From Appendix |A|, we see that the only relevant connection 
coefficient, T*^, is indeed equal to the corresponding Christoffel symbol. It therefore follows that Stpu 13 = -ft = E is 
a conserved quantity, again corresponding physically to the conservation of energy per unit mass. 

In contrast, consider the Killing vector X — d/d(j> of this same system. An inspection of the connection coefficients 
listed in Appendix [A] will show that these are not equal to the corresponding Christoffel symbols. It follows that 
s^pu^ = r 2 </>sin 2 9 is not conserved, despite the postulated spherical symmetry of the situation. This might lead one 
to believe that angular momentum would therefore not be conserved; in fact, it will be shown in the next section that 
there exists a conserved quantity whose asymptotic value is the conventional angular momentum. 

Let X 13 = u 13 in (0), giving 



//« q O«V- (8) 



In cases where the right-hand side of (||) vanishes, this yields another constant of the motion: the magnitude of 
the velocity. This will obviously be true for the geodesic equation, as the Levi-Civita connection is by definition 
compatible with the symmetric part of the fundamental tensor. In other situations, we are faced with the possibility 
that the magnitude of the velocity of a particle will not be conserved. Certain special cases arise when the chosen 
connection is not that of Levi-Civita, yet the right-hand side of (||) still vanishes. This is the case when considering 
a static, spherically symmetric background field. As the asymptotic properties of such a field in NGT force W M = 
|L3| , pO| , the right-hand side of (0) is seen to vanish, allowing the magnitude of the velocity to be conserved. However, 
it must be emphasized that this result is true only for a static, spherically symmetric field; non-Levi-Civita motion in 
a more generic NGT field would most likely leave the magnitude of the velocity unconserved. 

No discussion of conservation laws for particle motion would be complete without mentioning the matter-response 
equation, itself a conservation law. In NGT, the matter-response equation can be written |T^] 

d P K ~ Idxg^T^ = 0, (9) 

where 

t P x = \(g^ + g x ,T»»). (10) 

Given the nonsymmetric nature of both g^ v and T^ u , it is not obvious how ([h]) can be inverted to give T^ v in terms of 
t£. It will be sufficient for our needs to note that the relationship between t% and T^ v is independent of the properties 
of T^, and thus does not involve any acceleration terms. 

If we assume the simplest case of a monopole test-particle, we have 

T^cfx^O (11a) 



(x a - X a )T^ d 3 x = 0, (lib) 

and similarly for higher-order moments, the integration being carried out over a hypersurface of constant t, as per the 
definition of Papapetrou (see ||, p. 154). Here, X a is the position of the monopole, while x a is the position of the 
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observer. We take x° = X a = t and parametrize the motion using t. It is a simple matter to show that (|Tl|) imply 
similar multipole relationships for 

Using a method similar in spirit to that of Papapetrou (see ||, pp. 152-157), we can show that (^|) leads to 



where we define by — s Xu t^. The first term on the right-hand side of (|12j ) still contains T^ v ; however, given 
the relationship between T^ v and t^ v , we can assert that this term does not contain the acceleration of the monopole, 
and that this acceleration is contained only in the first term on the left-hand side. 

If T^ v is taken to be symmetric, we can show that r* 1 " = T* 1 "; this is the motivation for the definition of r^ v . It 
then follows that 

T ^d 3 :r = — — [ T m d 3 x, 
dt dt 



which causes the right-hand side of (12) to vanish. In such a scenario, we find that (JL2J) reduces to 



d f dX f) \ , 3 , dX^dX v „ 

— m— — +m\P } — — = 0. 

ds \ ds J X » vS ds ds 



where we have suggestively introduced 



Thus ( p^ ) reduces to the geodesic equation when rW vanishes. However, note that we have not required that g\p, v ] 
vanish: we are still working within the realm of NGT. 

Of course, T^ v ' — is a special case; indeed, (12) would seem to indicate that even a simple monopole particle 



couples directly to both T" 1 "' and g\n V ], and in NGT there is certainly no need to force such a contribution to vanish. 

The preceeding arguments would certainly seem to lend credence to the choice of the geodesic equation over the 
path equation as the correct description of motion in NGT; at the very least, ( p"2| ) is a strong argument for considering 
a geodesic equation with a supplementary coupling term on its right-hand side. However, the path equation has 
appeared sufficiently often in the NGT literature (see, for example, |ll| and jljj, and the references cited therein) to 
warrant the investigation of its physical predictions. In the remainder of this work, we will describe in parallel some 



of the properties of the geodesic and path equations that we have uncovered. In Section VIII below, we will return 
to the question outlined above: What happens if T^ 1 ^ ^= 0? 



IV. MOTION IN A STATIC, SPHERICALLY SYMMETRIC FIELD 

This section presents the equations and first integrals of ([[]) in the background geometry of the static, spherically 
symmetric solution described in Appendix |A|. These results will be used in the forthcoming sections to study radial 
and circular motions. We consider first the Levi-Civita connection, followed by the T-connection. 



We have already found in Section III that geodesic motion in a static, spherically symmetric field yields as constants 
of the motion the energy per unit mass E = yt, and the angular momentum per unit mass J = /?</> sin 2 8, as well as 
leaving k 2 = s a (3U a u^ unchanged. The 0-component of the geodesic equation, 

d 2 13' dOdr n n (d<f> 



dr 2 • Pdrdr ^^{Tt) = °' 

can be satisfied identically by choosing 6(jo) = tt/2 and 0(tq) = for some initial proper time To, corresponding to 
planar orbits. Using this, the r-component of the geodesic equation becomes 

d 2 r a' (dr\ 2 (3'J 2 jE 2 

dr 2 + 2a \dr) 2(3 2 a + 2a 7 2 ( ' 

Although this equation is identical in form to the corresponding equation in GR, the predicted motion is different 
owing to the different functional content of the metric functions. 
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Given these constants, the conservation of k 2 may be written 



2 a a (dry E 2 J 2 

Results similar to those listed above can also be derived for the path equation, the only difficulty arising from the 
slightly more complex nature of the connection coefficients. 



The i-component of the path equation was integrated in Section [II, and found to lead to a constant of the motion: 
E = "ft. This constant can be interpreted physically by examining its value asymptotically in the flat-space region, 
where it is found to correspond to the particle energy per unit mass. 



We have also discussed in Section III how, barring unusual circumstances, the path equation will in general not 
leave the magnitude of the velocity, k,~ unchanged. However, in certain special cases, of which the static, spherically 
symmetric field is one, a conserved k 2 can be extracted from the path equation. This is by virtue of the form of 
(||), whose rightmost member precludes any contribution from the antisymmetric piece of D a s^ u . It can be shown 
that, in the static, spherically symmetric solution to the NGT field equations, the only contribution coming from this 
covariant derivative is indeed antisymmetric on its indices. We therefore have a conserved k 2 : k 2 — s a f)U a u^ . 

As was the case for the geodesic equation, the ^-component of the path equation, 

d 2 9 pp' + ff drdO . Jd<f>\ 2 

can be satisfied identically by choosing 9(tq) — tt/2 and 9 (to) = 0, so that orbits will lie in a plane. Here, / is defined 
by f sin 9 = gwM- The (^-component of the path equation then leads to 



f± W + ff dr d4> = rfln(y/F+fW-/) = 
dr 2 (0 2 + f 2 ) dr dr dr 

where J is a constant having dimensions of a length. Solving for J, we have 



The physical significance of J is found by taking the limit r — ► oo. From (|A5|), we see that / — > exponentially, so 
that J behaves asymptotically as r 2 tfi, which we identify with the conventional angular momentum (per unit mass). 
We therefore conclude that it is J = \J ft 2 + f 2 <p and not 0(j> which represents the axial component of the angular 
momentum of the path equation. Near the origin, where / is not small with respect to 0, the value of J for a given 
cj) will be very different from the conventional angular momentum. 

Having determined the expression for the constant J, we can use this to give an expression for k 2 : 



^ = -s a ^ = a(^-) -^ + ^ (15) 



dr\ 2 E 2 J 2 (3 
dT~) ~ T + P 2 + f 2 

Use of these constants of the motion allows us to simplify the r-component of the path equation 

d 2 r a' fdr\ 2 i E 2 J 2 \3' {(3 2 - f 2 ) + 2/ /' 



dT 2 2a\dTj 2aj 2 2 a(0 2 + f 2 ) 2 



= 0. (16) 



We can combine path and geodesic motion by introducing the function i/ , ( r )j given by tp = ip g — 1/0 for the geodesic 
equation, and if> = ip p = 0/(0 2 + f 2 ) for the path equation. We then have 



dry , iE 2 jv =0 (17) 



instead of (Of) and (Rw, and 



instead of @ and @. 



dr 2 2a \dr J 2aj 2 2a 



K 2 = a(^-X~ — +J> (18) 
dr J 7 



G 



We can get a simple comparison of the geodesic and path equations by first expressing ( |17| ) in terms of r(t): 

d 2 r o^fdrV j (3 riln(V>) / d(j>\ 2 _ 
'dl 2 + 2^\dtj + 2a 2a dr \~db) 

The initial conditions for this differential equation are t = to, ro — r(to), ro = r(to), and </>o = 4>(to), where a dot 
represents differentiation with respect to t; note that we do not need to specify c£>q — 4>(to) as it does not appear in 
the differential equation. From the constancy of J and E, we have that 

PI Pa 

where we have adopted the notation that a subscript "0" on a quantity indicates that the quantity is to be evaluated 
using the initial data. It follows that the initial acceleration of the particle is 

d 2 r a r 2 7o | (3 <p 2 rilnQ) 

dt 2 q 2ao 2«o 2ao dr 

The rightmost term distinguishes the geodesic and path equations. It is always true that hi(ip g /ip p ) > 0, and for all 
r > M we have ln(M 2 ^ s ) < 0; therefore ln(M 2 ip g ) < ln(M 2 tp p ) < 0. In the large r limit, ln(-0 p ) ~ kx(i/) g ). If we 
assume that f(r) does not behave pathologically anywhere between r = M and r — > +oo, it then follows that 

d\n(ip g ) d\n(ip p ) 
dr ~ dr 1 

with dhi(ip p )/dr < 0. This allows us to conclude that initially, given identical initial conditions, a particle obeying 
the path equation would feel a greater inward radial acceleration than a particle obeying the geodesic equation. We 
could imagine a thought experiment where two particles, one obeying the path equation and the other the geodesic 
equation, are weakly deflected by some gravitational source. Given our previous analysis, we would expect that the 
first particle, in obeying the path equation, would be deflected through a larger angle than the second particle. In 
fact, although such a calculation would take us too far off course, an interesting method to determine experimentally 
whether an NGT test-particle obeys the path equation or the geodesic equation would be to set up just such a 
scattering experiment (see, for instance, plj , pp. I05-II9). Given certain values for the parameters, one might be 
able to determine a bound on the scattering angle with sufficient precision to distinguish between the geodesic and 
path equations. 

We will now use the results of this section to study the properties of radial and circular test-particle motion in a 
static, spherically-symmetric NGT background. 



V. RADIAL MOTION 



It is interesting to note that the geodesic and path equations differ only in their angular content: if J = 0, they 
predict the same motion. This can be seen in both ( [T^ ) and (18). It follows that a study of the radial motion of a 
particle is, in fact, a generic result, independent of which connection was used to study the motion. 

A particularly convenient form of the equation of motion is found by introducing reduced variables into Jl§|): 



dr x 2 



,ai- +7(1 + J4>) =K + V = E 2 , (19) 

where E — E/k, J — J/k and f = kt. K = ja(dr/df) 2 and V = 7(1 + J 2 ip) are the reduced kinetic and potential 
energies, respectively. 

From (p^), we can compute the acceleration of a radially infalling particle (J = 0): 

d f dr^\ 2 ^dr d 2 r E 2 d(aj) dr 1 da dr 



df \ df J df df 2 0:7 dr df a 2 dr df 
Dividing through by dr j df gives 
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d r 

d^ 2 



I (da E 2 d(aj) \ 
2a 2 \ dr j 2 dr J 



This expression is best evaluated in the ^-coordinate system (the ^-coordinate system is described in Appendix ^). 
It is found that 



dhj_ 
df 2 



2£(i/) 



2(asinh(ai^) + bsm(bis)) ~ 2 



where 



cosh(a^) — cos(6^) 

M 2 (l + s 2 ) 
(cosh(ai/) — cos(6t/)) 2 



(E 2 - e v ) 



(20) 



Since v < 0, we see that the particle is always attracted toward the larger negative values of v. 

The acceleration of a static particle is more easily found from (16) or (|l7|). Dividing this through by n 2 and setting 
dr/df = and J = gives 



d r 
df 2 



l V 2 



jE_ 

'2aj 2 



Again, the ^-coordinate system is more useful for evaluating this expression, giving 

E 2 



d 2 ^ 



2£(i/) 



(21) 



(22) 



To work out the turning points, we consider (|T^) in the ^-coordinate system; we find that there is only one turning 
point, at 



(23) 



where i>i is the initial radial position and is the initial radial velocity. In order to ensure that the particle velocity 
is real, we must have v < v t . However, since there is only one turning point, once the particle begins to head towards 
v — > — oo, it continues to do so, barring the application of some other force. It is therefore informative to compute 
the proper time required for the particle to reach some vf after starting from v%. This is found by inverting (|l^) in 
the ^-coordinate system: 



(M 2 (l + s 2 )) 



1/2 



dp 



(cosh(ai/) - cos{bv)){E 2 ~ e v ) x l 2 



We consider the case Vf < Vi < 0. Bounding E 2 — e v by E 2 — e Vi = £,(vi)vi 2 and cos(6^) by 1, we have that 

/M 2 (l + s 2 )\ 1/2 



f < 



V d 2 i(yi)i>i 2 



(coth(ai/{/2) — coth(af//2)). 



(24) 



This is a finite number for all Vi and Vf, with the exception Vi — > 0, where f — > +oo. We see that a particle can travel 
between any two non-zero v values in a finite amount of proper time. 

There is a further com plication associated with the physical interpretation of the results we have arrived at: it 
can been seen from (A4) that when v goes from v > vq to v < vq, the metric undergoes a signature change, with 
(3(v) going from being positive to being negative. At the point v — v$, the metric is necessarily degenerate. However 
we have found that a particle, regardless of whether it obeys the geodesic equation or the path equation, proceeds 
through the point v = vq and beyond without hindrance, despite the fact that at v = i/o, it would appear that the 
fundamental assumption of metric-non-degeneracy, so crucial to the postulates of the theory, is no longer valid. 

Although these results are reminiscent of the Schwarzschild singularity in GR (see 0], p. 663), it should be kept 
in mind that we are dealing here with a vacuum solution. It has yet to be established whether the Wyman solution 
is actually the final state of some generic stellar collapse problem (see p3| ) . 
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VI. CIRCULAR MOTION 



It is interesting to compare the predictions of the geodesic and path equations in the special case of circular motion 
as this displays the most striking difference between these two forms of motion: their ability or inability to support 
circular orbits for certain values of the parameters. 

In order to explore the existence of circular orbits, we minimize V with respect to r: 

_ dV_ = <h_ J2 rf (7^) 
dr dr dr 

Solving for J gives 

r2 _ d l/ dr r9C -N 

The location of the minimum of V is the radius of a circular orbit for a particular value of J. Typically in GR, the 
approach is to solve (|2^) for rasa function of J. Analysis of this equation shows that there is a minimum allowable 
value of the angular momentum, Jo. However, this approach is not practical in NGT, as the resulting equations are 
seemingly intractable. In turns out to be more profitable to turn the problem around and to view (Eq) as giving the 
value of J necessary to establish an orbit of a given radius. It is also useful to work in the ^-coordinate system, as this 



allows us to investigate the region near r = M. We can use (AS) to convert v values to conventional radial values. 

J is a one-to-many function of v. In particular, in the region v > z/ , J is a one-to-two function of v. There must 
therefore be a local extremum, which turns out to be a minimum. This local minimum will correspond to the smallest 
allowable angular momentum. We find this local minimum by solving 

dJ 2 d dj/dv 



dv dv d(^ip) I ' dv 



= (26) 



for v — f m , and then evaluating r m 2 — (3{v m ) and J 2 (y m ) = J m 2 . In GR, the resulting values (in conventional 
spherical coordinates) are r m — 6M and J m 2 — 12M 2 (see p2[ , p. 662) Unfortunately ( |26| ) cannot be solved in closed 
form when 7(f) and ij}{v) take on their NGT form ( |A2| ). However, in formulating the problem in this fashion, it 
becomes straightforward to solve ((2^) numerically for v rn . As an example, we take s = 0.9 and fxM — 0. Typically, 
a more reasonable choice might be to take \iM ~ 10~ 9 or so; however, since we are concerned with values of v near 
the origin, the problems generated in NGT by having fj,M = are not relevant here (see jl0| for a discussion of 
the problems associated with having \iM = 0), and taking \iM = simply represents an extreme case. The value 
s = 0.9 is chosen so as to emphasize the observed behaviour. Numerically, it is found that for the geodesic equation, 
v m w -0.40602, r m « 5.9930M, and J m 2 w 11.995M 2 , while for the path equation, v m w -0.40511, r m w 6.0040M, 
and J m 2 w 12.003M 2 . Generically, the geodesic equation yields smaller values than in GR, while the path equation 
yields values larger than the corresponding GR values. 

The advantage of using the ^-coordinate system is that the statements of the previous paragraphs are not tied 
down to some expansion, valid only in a limited region of spacetime; results calculated using the ^-coordinate system 
can be considered exact and valid for all v in the case of infinite-ranged NGT, or approximately true for finite-ranged 
NGT, in those regions where /ir < 1 and when /iM can be considered small. In particular, this is the case near the 
origin, where r ~ M <^ 1/fJi, and we may conclude that the path equation does not support circular orbits lying below 
r = 2M. On the other hand, circular orbits for the geodesic equation can be made to approach the origin arbitrarily 
closely by selecting successively larger values of the parameter s. 

Since the exact Wyman solution is not a solution to the field equations of massive NGT, one may question the 
wisdom of using it to investigate particle motion in NGT, not to mention the preceding treatment of circular orbits. 
Indeed, if applicable at all, the Wyman solution can only be used near the origin, where we would expect the range 
1 / n of the antisymmetric components of the field to have little effect. The study of circular orbits provides an excellent 
opportunity to show that using the Wyman solution gives qualitatively similar results to using an asymptotic solution 
to the field equations. In fact, if we repeat the preceding analysis of circular orbits in the path equation, but use 



(A5) instead of the Wyman solution, we find that the minimum circular orbit for s — 0.9 occurs at r m « 6.0033M, 
with J rn 2 m 12. 003-M 2 . The geodesic equation cannot be similarly verified, as its minimum circular orbits lie below 
r = 2M, where (|A5|) can no longer be trusted. 
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VII. BEHAVIOUR IN THE ASYMPTOTIC REGIONS OF SPACETIME 



By studying the geodesic and path equations in a static, spherically symmetric field, we were able to see in the 
previous two sections that there are differences in the resulting motion in the strong-field regime, in particular for 
the case of circular motion. In this brief section, we will demonstrate explicitly for the case of a static, spherically 
symmetric background field that the difference between the geodesic and path equations in the weak-field regime is 
of higher -ord er, and can safely be neglected. 



From (A5), we know that f(r) < r 2 for r > M. This behaviour is emphasized at larger r, as f(r) is damped 
exponentially Given this behaviour, we will consider the difference between the path equation and the geodesic 
equation by expanding in powers of £ — f(r)/r 2 ; we will show that this difference must be neglected in order to be 
consistent with our other approximations. 

Take ([!]) with C@ — and rewrite this as 

^ + {>^ = (0-r^)^ (27) 

by adding and subtracting a Christoffel symbol. 

Recall that the t- and ^-components had the same behaviour in both the geodesic and path equations. We need 
only consider the <fi- and r-components. Since the ^-component is immediately integrated, its expansion is trivial: 

d± ^ J 2 / _f 
dr r 4 \ r 4 

where J represents the value of the angular momentum in both the geodesic and path equations. Expanding the 
right-hand side of the radial component of (p7j) to lowest order in f(r)/r 2 , we arrive at 



d 2 r a' (dry j ' E 2 J 2 fff 3/ 2 \ J 



dr 2 2a \dr J 2aj 2 or 3 \ r e r 7 J a 



(28) 



The left-hand side of this is nothing more than the geodesic equation. The right-hand side is of order e 2 ; approximating 
f(r) and f'(r) by their maximum values of f(r) ~ sM 2 /3 and f'(r) ~ — s/iM 2 /3, we could conclude that, for large r, 
the difference between the path equation and the geodesic equation is of order 1/r 6 . Adding to this the exponentially- 



damped behaviour of /(r), we could conclude that for any r > 2M where (A5) can be trusted, the geodesic and 
path equations will yield identical results. Much stronger than this, however, is the fact that the corrections to 7(f) 
and a(r) in QAjl) are themselves of order e 2 (see (4.6) in [Of, and the discussion preceding it). Thus, if we accept 



(A5) as valid approximations to some exact solution of the NGT field equations, in order to be consistent with this 
approximation, we must drop the right-hand side of (p8|). 

This motivates the use of the geodesic equation in weak-field regions of spacetime. 



VIII. MODIFYING GEODESIC MOTION 



In the previous section, it was shown that path motion converges very rapidly to geodesic motion. However, by its 
very nature geodesic motion does not couple the particle directly to the antisymmetric components of the NGT field. 
Moreover, we saw at the end of Section III that even a monopole test-particle must, in the most general case, couple 
to a^ v . In this section, we present three possible couplings. For simplicity, these are linear in the particle velocity. 

Since we are dealing in this section only with the weak-field regions of spacetime, we will consider modifying only 
the geodesic equation; from the previous section, we know that the path equation will yield similar results in these 
regions. 

The couplings will be derived from a scalar Lagrangian L, which we take to have the form L = (l/2)A fl UL IJ ', where 
A is a covector independent of the particle velocity. The Euler-Lagrange equation of particle mechanics allows us to 
conclude that the contribution of such a Lagrangian to the equation of motion will be of the form 



where f [a ^ = d [a A^. 

Let L = L\ + L 2 + £3, where 
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1 „ Ai e wl/A 

Li = -Ai*F''u r; = — j— F[ lluX ]S V aU ,y (29a) 

£a = ^A 25 [^lF [Ml/A]U A (29b) 

L 3 = ^M=- aA/i „P (29c) 



2 V=5 

The components of the field-strength tensor F^x] arc defined by 

F[nu\] = d{\g^u] = ^(dxa^u + <9 M a„ A + d„a x ^). (30) 

The constants Aj (i — 1,2,3) couple the test-particle to the NGT skew field, and have dimensions of a length. The 
symbol 71 is the fully antisymmetric Levi-Civita tensor density, defined by 

T if iivXq is an even permutation of 0123, 
e ^u\i) _ ^ _^ jf ^ v y^ j s an Q(j(j permutation of 0123, 

otherwise. 

It will be found below that only L% will generate any interaction at all in a static, spherically symmetric field. This 
does not, of course, exclude the possibility that L 2 and L 3 could more accurately reflect the desired coupling between 
the test-particle and the antisymmetric components of the NGT field. It does, however, make their contribution 
trivial in our particular case of interest. Although for the sake of completeness we will derive the additional terms 
in the equation of motion for the case when all three interactions are present, we will in the end ignore Li and L3, 
knowing that they will in no way affect our results. 

Using the Euler-Lagrange equation, we find that 

f[ap\ = X ld[ a y r—^ F[<ru\]9n]r)J + ^2 d [cx \ 9 F [wn]]) + X sd[ a I ^— S^j , (31) 

such that the equation of motion is now written 



dT 



+ CH^ = ^ a f[a^ (32) 



where k 2 = s^u^u 1 ' is the magnitude of the veloci ty. 

As we will only be interested in the behaviour of (B2h in a static, spherically symmetric field, it will be of no interest 



to re-derive the conservation laws discussed in Section III. Suffice it to say that the techniques of that section can 



easily be generalized, and the first integrals are found to be slightly different. It will be simpler for us to derive the 
constants of the motion "by hand" below. It is straightforward to show that k 2 is conserved. The appearance of 
the magnitude of the velocity in the equation of motion in this fashion has an important implication: for a massless 
particle (such as a photon), k 2 — and the right-hand side of (|32|) is seen to vanish. Therefore, in this scheme, 
massless particles will not couple to the antisymmetric field. 

The case n 2 = corresponds to pure geodesic motion, which was discussed in previous sections. We will therefore 
concentrate on the case k 2 ^ 0. Since k 2 is conserved, we can scale the proper time r in such a way that k — 1; 
henceforth, we assume that this has been done. 

In the static, spherically symmetric field, the skew field-strength tensor -Fj^Al nas onr y one independent, non-zero 
component, 

F[e<t>r] = ^9 r ae4, = ^/' sin61 - 

On the other hand, because of the chosen form of the fundamental tensor, d„g[ Al/ l can be shown to vanish. From (|3 
it follows that the tensor jy aG \ also has one independent component: 



Ai7/' 



» \^7(r 4 + / 2 )/ 

For convenience, we will write Ai = A, since the A2 and A3 terms in ( |3l| ) make no contribution to the equation of 
motion. 
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Using (|32[), we find that the equations of motion for a particle in the field of a static, spherically symmetric source 
are given by 



= 

= 
= 



dH 
dr 2 

d 2 ^ 
dr 2 

fl 

d^ 2 
d?0 
d^ 2 



7' dt dr A dr d 
7 dr dr 7 dr dr 



if 



a 



drY 
2a \dr J 



\/aj(r 4 



dr 



P)J 

d6 
a \ dr 



1 



dtY 
2a \dr) 



A dt d 
a dr dr 



if 



V^lir 4 + P) 



2 d6 dr 
r dr dr 



sm f cos f 



dr J 



<d0 



2 dch dr cos t 

- — h 2- : . 

r dr dr sin 9 dr dr 



(33a) 

(33b) 

(33c) 
(33d) 



We can satisfy (33c) identically by letting 9(tq) = ir/2 and 9{tq) = for some initial proper time To. It follows that 
orbits lie in a plane. Meanwhile, we conclude from ( |33d| ) that J — r 2 <j) is a constant of the motion. This might have 
been predicted from the fact that s^ a /[ ajU ]U M , the right-hand side of (|32|), does not have a 0-component, and as was 
seen in Section III, J is conserved for geodesic motion. 

Contrary to the situation for geodesic and path motion, s t pu@ = ^yi is not conserved. In fact, we see from (33a) 
that the actual first integral of the motion is 



dt 

= 1 dr 



A7/' 



vW 4 + P) 



(34) 



That this represents the energy per unit mass can be seen by taking the large r limit and using ( |A5| ) . 
With these results, we can rewrite ( 33b ) as 







fr 
dr 2 



dr 
2a \dr 




A7/' 



V«7(r 4 + P) 



(35) 



There are two calculations to perform. First, we calculate the corrections to the Newtonian gravitational force 
coming from the supplementary forcing term. Secondly, we will determine the equation for the orbit of a particle 
acted upon by this extra force. In order to extract a useful result from this, we will make certain simplifying 
assumptions. In particular, we will assume that f(r)/r 2 <C 1 and \f'(r)/r 2 <C 1. If s <C 1 and /J.M <C 1, this 
corresponds to regions of spacetime where Mjr 1. 

Consider first the Newtonian gravitational force. We rewrite (^5|) in terms of r(t): 




if 



dr 



VaM^+p) 



Here, Jn = r 2 d<f>/dt is the Newtonian value of the angular momentum per unit mass. We recognize the left-hand side 
as the usual GR contributions to the equation of motion. 
To our order of approximation, 

XsM 2 fi 2 e~^ r (l + nr) 
3~E V 2 ' 




Therefore, the radial equation of motion may be written 



cPr 
dt 2 



Jn 



M 



XsM 2 fi 2 e-' tr (l + fir) 
3~E r 2 



(36) 



where we have assumed that the particle is moving slowly, so that dr / dt < 1. If A < 0, this yields a repulsive Yukawa 
force in ([36]), while A > renders this force attractive. 

The Yukawa behaviour generated by this modification of geodesic motion has interesting phenomenological impli- 
cations. Take, for instance, the case of very large r. We find in this case that the right-hand side of (p§) reduces to 
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—M/r 2 ; that is. the usual Newtonian law of gravitation. On the other hand, for r -C 1/fJ, but still with r 3> M, the 
last term in (13q) becomes 



1 



3E 



M 



Again, this is the usual Newtonian law of gravitation, but with a renormalized gravitational constant. Therefore, in 
the two weak-field regions that we have just considered, we find that Newtonian gravitation applies. However, the 
actual strength of the force is different in each region. This can be used to fit the galaxy rotation curves without 
resorting to postulating dark-matter [^jj . 

For completeness, we will now derive the equation for the orbit of a particle in this scheme. From the normalization 
condition s^u^u" = 1 and using 8 = 7r/2 and 6 = 0, we arrive at 



du 



1 



E 



At/' 



i 



i 



V«7(r 4 + P) J 

Differentiating this with respect to <fi gives the equation describing the orbit of the particle: 

At/' V 



d(f> 2 



ajJ 2 



E- 



V«7(r 4 + f 2 ) 



1 



(37) 



Thus far, the treatment has been entirely general. We can specialize this to t he case where f(r)/r 2 <C 1 and 
Xf'(r)/r 2 <C 1. Taking the Schwarzschild forms 7 = 1/a = 1 — 2M/r and using (A5), we arrive at the orbit equation: 



I2 



M , 

T 2{] 



EXsM/j, 2 
12 



(1 + \ir)e~ 



MIu z 



(38) 



This equation is similar to the orbit equation from GR (see p5| , p. 186), the only difference being the factor multiplying 
the angular momentum term on the right-hand side. 

A comment should be made about the role played by the quantity A in these equations. In a sense, A represents 
the strength of the coupling between the test-particle and the NGT field; it could be thought of as a charge, although 
unlike the electromagnetic charge, it does not follow from any sort of flux integral. The question immediately arises 
as to whether A is a property of the test-particle or some universal constant. If the former, then we must deal with 
the violations of the weak equivalence principle that this material-dependence will cause (see |2f|, p. 13). However, 
since there is no compelling evidence to suggest that the weak equivalence principle is actually violated in reality [p7| , 
it would seem more sensible to choose the latter route and take A to be some material- independent constant. 



IX. CONCLUSIONS 



We have described observational differences between the geodesic and path equations, in the hope that experiments 
may differentiate between the two. However, despite the differences between these two equations of motion, we 
have found that under certain circumstances, for example the case of radial motion, they predict identical results. 
Furthermore, we have shown that in the weak-field limit, both of these equations lead to virtually identical equations 
of motion, owing to the fact that NGT is constructed so as to recover GR in the weak-field limit. 

Many of our results are implicitly linked to the case of a static, spherically symmetric field. This is a question of 
practicality: as it stands, this is the only available solution to the field equations. However, we would expect that 
the qualitative effects that we have described would be generic. In particular, as pertains to constants of the motion 
and the like, the geodesic equation of NGT behaves very much like its counterpart in GR. This is to be expected: the 
geodesic equation of NGT is constructed so as to couple the particle only indirectly to the antisymmetric field. On 
the other hand, it was found that the path equation had different constants of the motion. In particular, we found 
that the conventional angular momentum, J = 13(f) sin 2 6, was not a first integral of the motion. Rather, a different 
first integral was derived, and was found to be a sort of "generalized" angular momentum. Although we have no proof 
to support this, the work we have performed thus far would lead us to expect that this generalization of constants 
of the motion would be a generic feature, owing to the right-hand side of (^): whenever we have come across a case 
where (^) did not yield an explicit constant of the motion, the resulting equation could nonetheless be integrated to 
yield a first integral. 



13 



Particularly interesting behaviour was found by studying radial motion in a static, spherically symmetric field. 
Firstly, it was found that an infalling particle need not stop at the origin r — 0. This was seen using the ^-coordinate 
system, where it was found that a particle can pass through vq and proceed toward v < i/g. Furthermore, the particle 
could reach any radius Vf from any non-zero radius Vi in a finite proper time. An upper bound was placed on the 
proper time taken to travel this distance. It is unclear at this time exactly what is implied by this result: the Wyman 
solution to the infinite-ranged NGT field equations is a vacuum solution, devoid of any matter contributions. As 
the gravitational collapse problem in NGT remains unsolved, it is not obvious that the vacuum Wyman solution, as 
opposed to, say, the Schwarzschild solution, is in fact the final state of a stellar collapse. Regardless, the comparison 
with the corresponding behaviour in GR is striking, and should be kept in mind. 

The major difference between the geodesic and path equations was found by studying circular motion. Here, it was 
shown that circular orbits for the geodesic equation in NGT extend to lower radii than in GR, while the last circular 
orbit for the path equation lies outside its GR counterpart. This behaviour was found to be strengthened for larger 
values of the parameter s appearing in the static, spherically symmetric field. 

After demonstrating that the geodesic and path equations have similar weak-field limits, we introduced three 
interactions which could serve as alternate methods for coupling a test-particle to the antisymmetric components of 
the NGT field. Of these three, only one was found to generate any interaction at all in a static, spherically symmetric 
field. The correction to the Newtonian gravitational force acting on the particle was worked out in the weak- field limit, 
and was found to be a Yukawa force. The correction to the orbit of a particle was also calculated; the lowest-order 
NGT correction was shown to be a l/r 3 term. 
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APPENDIX A: THE STATIC, SPHERICALLY SYMMETRIC SOLUTION IN NGT 



We list here the salient features and results of the static, spherically symmetric solution to the new NGT field 



t, x 1 



r, x 



6> and x 



where r, Q, and d> are the usual 



equations. The coordinate system used is x° 
spherical coordinates. 

The NGT field equations contain a parameter 1 j \i which is identified with the range of the skew components of the 
fundamental tensor. The static, spherically symmetric solution introduces two constants of integration: M and s. M 
plays the same role as the Schwarzschild mass in GR, while s is a dimensionless constant. Beyond being a measure of 
the strength of the antisymmetric contributions to the gravitational field, the physical meaning of s is still unclear. 
To be specific, we will assume that M <C 1/m- 

The fundamental tensor is written 



9tw 



7 (r) 


w(r) 








—w(r) 


—a(r) 














-/?(r) 


f(r) sin 








—f(r) sin 6* 


-^sin 2 



(Al) 



In M] and |2Q], it is shown that the only solution which yields an asymptotically Minkowskian space has w(r) = 0. 

In the case of infinite-ranged NGT, corresponding to setting fj, = in the field equations, an exact solution to the 
field equations exists: 



7(r) 
a(r) 
/9(r) 

m 



M 2 e~ u (l + s 2 ) 
(cosh(a^) — cos(&^)) 2 

„2 



2Af 2 (sinh(ai/) sin(&^) + s(l — cosh(ais) cos(bu))) 
e"(cosh(aj/) — cos(fe^)) 2 



(A2a) 
(A2b) 
(A2c) 
(A2d) 
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where 



and 



The function u{r) is determined from the relation 
e u {cosh{av) - cos{bv)) 2 = 



cosh(a^) cos(fe^) — 1 + s sinh(ai^) sin(6f). 



(A3) 



This exact solution is referred to as the Wyman solution |28||29|], and holds for all s. 

There are two coordinate syste ms o f interest in the Wyman solution. The first is the conventional set of spherical 
coordinates used above. However, ( |A3| ) can be viewed as defining a transformation from a different coordinate system, 
the ^-coordinates, to conventional spherical coordinates pCj ], The ^-coordinates are interesting, as (AS) need not be 
inverted; the coordinate system is x° = t, x = v, x 2 = 6 and x 3 = </>, where 



a{v) 



M 2 e-»(\ + s 2 ) 
(cosh(a^) — cos(fe^)) 2 



2Af 2 (cosh(ai/) cos(fr^) — 1 + s sinh(ai') sin(6i/)) 
P{y > = e-(cosh(ai/)-cos(6^)) 2 ' ( A ' 



7 (r) s 


. 1 2M 


a(r) s 






3 



and with 7(f) and f(y) given as above. Since (A3) gives r as a many-to-one function of ^, a particular branch of the 
solution must be picked. This selection is done by picking the branch that will yield the positive-mass Schwarzschild 
solution as a limit. This branch begins at v = and extends toward negative v. In such a coordinate system, the 
asymptotic, weak-field region is at v — 0, while the "origin" occurs at vq defined by /3(fo) = 0. The particular value 
of vq depends on the value of s; for s — 1, we find numerically that vq ~ —5.1667. 

When f{r)/r 2 <C 1, a solution to the linearized field equations is approximated by pc| ] 

(A5a) 
(A5b) 

■fir); (A5c) 

the coordinates are the conventional spherical coordinates. In [ jl3| , a similar result is derived. The two forms can be 
seen to be asymptotically equivalent, notwithstanding a change in notation, by taking fiM <C 1 in (2.12) of Jl3| and 
taking fir 1 in the above result. 

For both conventional spherical coordinates and the ^-coordinates, the Christoffel symbols are given by: 

{/j = i, (A6a) 

{ X u) = Ya (A6b) 
I x i = — (A6c) 
{u)=^ 2 ^e} = -^ 1 (A6d) 

1%} = l\) = § (A6e) 
Q =- sin0 cos 9 (A6f) 

«>-££ (A6g) 

(A6h) 
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where x represents either r or v. A prime denotes differentiation with respect to x. 

The T-connection coefficients, are determined by solving the NGT compatibility condition: 



(A7) 



where is determined from 

d P g [ap] = -\g (pa) W p . 

When w(r) — 0, it can be shown that W M = 0. The right-hand side of ( |A7j ) therefore vanishes. For both conventional 
spherical coordinates and the ^-coordinates, the nonvanishing components of the T-connection are: 



pt / 

1 (tx^ — ~ 



(tt 



pa; 

(xx 



-px 



(x6 



r 
r 

r 

n 



(OA 



2_ 

= L 

2a 

- ^1 
_ 2a 

= sin 2 ^r? flfl , = - 



(/3'(/3 2 -/ 2 ) + 2//'/3)sin 2 



2a(/3 2 + p) 

0(3' + ff 



r (x<t>) 2(^2 + /2 ) 

= — sin 9 cos 6* 
sin0 



cos 9 

(/'(/? 2 -Z 2 )- 2/3/3'/) sin t 



2a(/3 2 + J 2 ) 
■sm J 0r? on = 



{pf-pf)wH 



2(/3 2 + / 2 ) ■ 

Again, x denotes either r or v, and a prime denotes differentiation with respect to x. 



(A8a) 

(A8b) 

(A8c) 

(A8d) 

(A8e) 
(A8f) 
(A8g) 

(A8h) 
(A8i) 
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